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Abstract. The nonlinear Schrodinger (NLS) equation is a fundamental model 
for the nonlinear propagation of light pulses in optical fibers. We consider an 
integrable generalization of the NLS equation which was first derived by means 
of bi-Hamihonian methods in [A. S. Fokas, Phys. D 87 (1995), 145-150]. The 
purpose of the present paper is threefold: (a) We show how this generalized NLS 
equation arises as a model for nonlinear pulse propagation in monomode optical 
fibers when certain higher-order nonlinear efi^ects are taken into account; (b) We 
show that the equation is equivalent, up to a simple change of variables, to the 
first negative member of the integrable hierarchy associated with the derivative 
nonlinear Schrodinger equation; (c) We analyze traveling-wave solutions. 
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1. Introduction 
The nonlinear Schrodinger (NLS) equation 
(1.1) iut + ^Uxx + p\u\^u = 0, X G M, t > 0, 

where 7 and p are real parameters and n(x, t) is a complex-valued function, is a 
fundamental model for the nonlinear propagation of light pulses in optical fibers 



[T]. Equation (1.1) is a completely integrable system and can be linearized by the 



inverse scattering transform formalism pi. In the context of fiber optics, the soliton 



solutions of ( 1.1 ) are of particular interest: because of their inherent stability, solitons 
can travel unchanged over long distances, making them ideal carriers of information 
in optical transmission systems. In this paper we consider a generalization of the NLS 
equation, which also admits soliton solutions, and which, as we shall argue, models 
the propagation of nonlinear light pulses in optical fibers when certain higher-order 
nonlinear effects are taken into account. 
The equation 

(1.2) iut — I'Utx + lUxx + p\u\^{u + ii^Ux) = 0, X € M, t > 0, 

where z/, 7, p are real parameters and u{x, t) is a complex-valued function, was derived 
in [3j as an equation with two distinct, but compatible, Hamiltonian formulations. 
Just like the bi-Hamiltonian structure of the well-known Korteweg-de Vries equation 
can be perturbed to yield the integrable so-called Camassa-Holm equation |4j, the 
same mathematical trick applied to the two Hamiltonian operators associated with 



the NLS equation yields equation (1.2) — note that \1.2\ reduces to (1.1 ) when 1/ = 0. 



Equation (1.2) admits a Lax pair formulation and the initial value problem can be 



solved by means of the inverse scattering transform [5] . 
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For large classes of dispersive nonlinear PDE's, the NLS equation arises asymptot- 
ically as the equation satisfied at first approximation by the amplitude j4 of a slowly 
modulated monochromatic wave [6]. It was demonstrated in [3] that an application 
of asymptotic techniques to a KdV type equation gives rise to an equation of the 



form (1.2) provided that one retains terms of the next asymptotic order beyond those 



included in the case of the NLS equation. In this sense equation (1.2) is a general- 



ization of the NLS equation. However, to the best of our knowledge a more direct 



physical derivation of (1.2) has not yet been presented. 



The purpose of the present paper is threefold: (a) We show how equation (1.2) 
arises as a model for nonlinear wave propagation in monomode optical fibers; (b) 



We show that (1.2) is equivalent, up to a simple change of variables, to the first 



negative member of the integrable hierarchy associated with the derivative nonlinear 
Schrodinger (DNLS) equation; (c) We analyze traveling-wave solutions. 
Let us comment on (a), (b), and (c). 

(a) Physical derivation. The propagation of nonlinear pulses in optical fibers is 
described to first order by the NLS equation (1.1). However, when considering very 



short input pulses it is necessary to include additional terms accounting for higher- 
order effects such as self-steepening, Raman scattering, and third-order dispersion 
[JIIH]- While the rmore general esulting equation is not integrable, there exist exactly 
solvable approximations which capture various aspects of the higher-order phenom- 
ena. Most notably, the DNLS, Hirota, and Sasa-Satsuma equations are all integrable 
reductions of the more general equation and are used as models when studying wave 
propagation for ultrashort input pulses O [71 |T0]. Although each of these reduc- 
tions ignores some of the important nonlinear effects, their solutions can be used as a 
first approximation when performing a more complete analysis involving perturbation 
methods. Moreover, the possibility of presenting exact solutions to these models make 
them valuable tools for understanding the physical influence of the various terms. 

After recalling certain aspects of the standard derivation of the NLS equation in 
nonlinear fiber optics, we indicate in Section [2] how equation (1.2) appears when 
taking into account terms that are normally ignored. 

(b) Relation to DNLS. The derivative nonlinear Schrodinger equation 



:i.3) 



iqt + Qxx + cri{\q\'^q)x = 0, a = ±l, 



and gauge transformed versions thereof, have several applications in plasma physics 
and nonlinear fiber optics |1H [TO] . It was shown in [12^ that (1.3) admits a Lax 
pair formulation and can be solved by means of inverse scattering techniques. Be- 



ing integrable, equation (1.3) admits an infinite number of conservation laws and a 



bi-Hamiltonian formulation. Moreover, associated with equation (1.3) is an infinite 



hierarchy of equations generated by the bi-Hamiltonian structure, equation (1.3) be- 



ing the second positive member. In Subsection 3.2 we demonstrate that, up to a 



simple change of variables, equation (1.2) is nothing but the first negative member of 
this hierarchy given by 



(1.4) 



Utx 



a 



±1. 



Since the parameter a in (1.4) can be removed by the replacement {t,x) {at, ax) 



this shows that the parameters z/, 7, p in equation (1.2) can all be set to 1 by a 



change of variables. In particular, in contrast to the case of the NLS equation which 
comes in a focusing as well as in a defocusing version depending on the values of the 
parameters, and solitons only exist in the the focusing regime, all versions of equation 
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(1.2) are mathematically equivalent up to a change of variables. This observation is 



consistent with the fact that solitons for (1.2) exist for all values of z^, 7, p |5j. 



Our investigation of the relationship between ( 1.2 ) and ( 1.3 ) was motivated by the 



observation that if Ux is identified with q, then the Lax pairs associated with equations 



(1.2) and (1.3) have identical x-parts [51. Let us point out that a similar link exists 
between the KdV and Camassa-Holm equations: The Camassa-Holm equation is 



derived mathematically from the KdV equation in an analogous manner that (1.2) 
is derived from the NLS equation. Moreover, it is known that the Camassa-Holm 
equation is related by a (nonlocal) change of variables to the first negative member 



of the KdV hierarchy [13]. Since the gauge transformation relating (1.2) to (1.3) is a 
local change of variables, the link in the present case is more direct. 



(c) Traveling waves. In Section |4] we consider traveling-wave solutions of (1.2). 
For simplicity, we consider the simpler, but equivalent, equation (1.4). If we assume 



that the solution u{x,t) of (1.4) has the special form 

(1.5) u{x, t) = ip{x - ct)e'^kx-nt+e{x-ct)) ^ 

where k, O, c are real parameters and and 6 are real-valued functions, the partial 
differential equation for u reduces to a first-order ODE for ip. A straightforward anal- 
ysis of this equation yields the qualitative structure of large classes of traveling-wave 
solutions. For certain values of the parameters the equation for ip can be integrated 
explicitly. In particular, we recover the one-soliton solutions earlier found by means 
of the inverse scattering approach [5]. We also comment on the occurence of more 
exotic traveling waves such as peakons. A peakon is a solution whose profile has a 
peak at its crest, see Figure [3j The existence of peaked traveling- wave solutions is 
a well-known property of the Camassa-Holm equation [14j. Since (1.2) is related to 



the NLS equation by a procedure analogous to that which gives the Camassa-Holm 



equation from KdV, it is natural to investigate whether (1.2) also exhibits weak solu- 
tions of this kind. Our discussion in Section |4] suggests that (1.2) admits no peaked 
traveling waves which are weak solutions in any reasonable sense. 

2. Physical derivation 

Consider the propagation of an optical pulse in a monomode fiber aligned in the z- 
direction with a frequency-dependent dielectric constant. We first follow the standard 
derivation of the NLS equation (see e.g. [1]), before we indicate how the additional 
higher-order effects particular to equation (1.2) arise. Maxwell's equations lead to 
the basic equation 



(2.1) 



V X V X E 



.E. 



where E(r, t) and P(r, t) are the electric field and the induced polarization evaluated 
at time t at the spatial point r, respectively, (j,q is the vacuum permeability, c is 
the speed of light in vacuum, and subscripts denote partial derivatives. Assuming 
a local medium response and including only third-order nonlinear effects, the total 
polarization P(r, t) = ^ li'^-, t) + 'Pml^'"^, t) can be split into a linear part 



X 



(1) 



{t-t')E{r,t')dt', 



and a nonlinear part 
ft ft 
'Pnl = eo / dti 



f dti I dt2 [ dt3x'-^\t-h,t-t2,t-t3)[E{r,ti),E{^,t2),E{r,ts)]dt', 
J —00 J —00 J —00 
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where cq is the vacuum permittivity and x^"'^ is the jth-order susceptibihty. If we 
moreover assume that the nonhnear response of the fiber is instantaneous, we find 

(2.2) P,vL(r,t) = eoX^'^[E(r,t),E(r,t),E(r,t)], 

where x^'^'^ is a time-independent map Unear in its three arguments. Since we consider 
a quasi-monochromatic pulse which maintains its polarization along the fiber length, 
a scalar approach is appropriate. If the pulse spectrum is centered at ujq, we write 

E(r,t) = |(i?(r,t)e-^'^»* + c.c.), 

(2.3) Pz.(r,t) = |(Pi(r,t)e— o* + c.c.), 

PiVL(r, t) = I {PNL{r, t)e~'^^' + c.c.) , 

where x is a unit vector pointing in the direction of the polarization, c.c. stands for 
the complex conjugate, and P^, P/vL are complex-valued slowly-varying functions. 
Substituting (2.3) into (2.2) and ignoring terms of frequency 3u;oj we deduce that 
Pnl = ^o^nlE, where e^L = ^x'xxxx\E{r,t)\'^ is the nonlinear contribution to the 

(3) 

dielectric constant responsible for the Kerr effect ( xxxx x denotes the xxxx-component 
of T^(3)). Taking the Fourier transform of equation (2.1 ) while making the approximate 
assumption that e^i is a constant, we obtain 

(2.4) V'^E{r, uj-ojo) + e{uj)klE{r, w - wq) = 0, 
where ko = co/c, the Fourier transform E is defined by 

E{r,Lo-uJo)= I E{r,t)e'^'^''^"^^dt, 
Jr 

and the frequency-dependent dielectric constant 

contains the Fourier transform (cj) of the xx-component of x''^''- 
We seek a solution to (2.4) of the form 

^(r, LO-uJo)= F{x, y)A{z, lo - ooo)e'^^' , 

where A{z, lo) is a slowly varying function of z and /3o is the wave number accounting 
for the fast oscillations. Substituting this into (2.4), we find 

(2.5) F^^ + Fyy+{e{uj)kl-~P{ujf^F = ^, 

(2.6) A,, + 2i/3oA, + (/3(u;)2 - /3g) i = 0, 

where /3(a;) is a function independent of r. At this point we have deviated from the 
standard derivation in which the A^z term in equation (2.6) is ignored since A{z,t) 
is assumed to describe the slowly varying envelope of the pulse. 

The derivation proceeds by solving (2.6 ) using perturbation theory and then apply- 
ing the inverse Fourier transform in order to return to the time domain. The result 
is (see [1] for details in the case when A^z is ignored) 

(2.7) - ^^Azz + A, + PiAt + '-^Att + f ^ = ip\A\^A 

where the terms proportional to a and p include the effects of fiber loss and non- 
linearity, respectively. The coefficients Po,Pi,P2 are determined by the expansion of 
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(3{uj) = n(uj)uj/c, where n{uj) is the refractive index, around the carrier frequency ujq 
according to 

/?M = ^0 + Pi{u^ -uo) + ^{uj- ujof + ^(a; - c^o)' + • • • • 

In the frame of reference travehng with the pulse at the group velocity Vg = l//3i, 
equation (2.7) becomes 

(2.8) iA, + :^A,, - -^A,T + iAtt + = -p\A\^A, 

where 7 = 2/3^ ~ ^ have introduced the new variable T hy T = t — Piz. 

As compared with the standard derivation, we see that our retaining of Azz in 
equation (2.6) has given rise to two additional terms in (2.8) involving A^z and Azt- 
If these are ignored and the dissipation is set to zero (i.e. a = 0), equation (2.8) 
reduces to the NLS equation. 

The terms proportional to Azz and Azt sue usually ignored due to the inequalities 
(see |15j) 

\A,z\ <C ^o\Az\, \Azt\ < PoVg\Az\, 

assumed to be valid if ^ is a slowly varying envelope. However, in the (femtosecond) 
regime of very short pulses, the pulse envelope A may contain only a few optical cycles. 
Hence these inequalities are expected to be less strong in this range and it motivates 
us to not throw away the terms involving Azz and Azt- On the other hand, for very 
short pulses the approximation (2.2) also has to be corrected in order to incorporate 
higher-order nonlinear effects. A more careful analysis [8j taking additional effects 
into account shows that the right-hand side of equation (|2.8|) should be replaced by 



'-^Attt - pA\A\^ - is {A\A\^)^ - irA {\A\'')^ , 



where /Js governs the effect of third-order dispersion, s accounts for the so-called self- 
steepening effect, and r is in general a complex parameter arising from the retarded 
nonlinear reponse of the medium; in particular, the imaginary part of r describes 
the retarded Raman effect in which a photon through scattering pass on part of its 
energy to a vibrational mode of a molecule in the medium. The Raman effect is a 
nonlinear dissipative effect which causes a frequency downshift of the pulse [U |8] . 
While our final equation 

2.9 lAz + -^r^Azz - -5 — AzT + iAtt + ^A- -^Attt 

2po PoVg 2 6 

= -pA\A\^ - is {A\A\^)^ - irA {\A\^)^ 

is rather complicated, there exist some special choices of the parameters for which the 
equation is integrable and hence exactly solvable by means of the inverse scattering 
formalism. First of all we need to set the linear and nonlinear dissipative terms to 
zero (i.e. a = and Imr = 0) in order to get an energy-conserving system. Then 
several integrable reductions exist. 

• When the terms proportional to Azz, Azt, Attt, and A (| Ap)^ are neglected, 
equation (2.9) reduces to 

(2.10) iAz + -/Att = -pA\A\^ - is {A\A\^)^ . 

This equation is integrable and admits soliton solutions [16]. In fact, a gauge 
transformation converts it into the DNLS equation (1.3). 



6 



JONATAN LENELLS 



(2.11) 



More generally, when the terms proportional to Azz and Azt are ignored, we 
obtain 

' -^'^'2 is{A\A\^)^-iTA{\A\')^. 



iAz + 7^TT 



6 



-A 



TTT 



-pA\A\' 



There are two known [T7] integrable subcases of this equation apart from 
(2.10); both occur when the parameters satisfy P^p = —2s'y. When this 



(2.12) 



condition is fulfilled, a gauge transformation brings ( 2.11[ ) to the form [I 



iA,, 



6 



'-Attt = -is\AYAT - i{s + t)A [\AY) 



T • 



If the ratios of the coefficients — s, and s+r in (|2.12|) are related as 1 : 6 : 



we obtain Hirota's equation |19j 

Az + Attt = -Q\A\'^At, 
while if they are related as 1 : 6 : 3 we find the Sasa-Satsuma equation [2U] 



Az + Attt = -Q\A\^At - {\A\'^) 



T • 



(2.13) 



When the terms proportional to Azz and Attt are ignored and s + r = 0, 

3S 

AzT + ^Att 



equation (2.9) becomes 
1 



iAz 



PoVg 



-P\A\ 



^ ' A + i-AT 
P 



which is equation (1.2) if 



1 



We conclude that these four exactly solvable models arise as reductions of the full 



nonlinear equation (2.9) when different choices are made as to which higher-order 



effects to include. Although none of them can be expected to exhibit all the structure 



tied into (2.9) (e.g. they all ignore the effects of fiber loss), they do provide important 



first approximations to the full theory. In this sense equation (1.2) can be viewed as 
a model for light pulses in single-mode optical fibers. 

3. Relation to DNLS 



It was found in |5j that the x-part of the Lax pair of (1.2) is simply related to the 
a;-part of the Lax pair of the derivative nonlinear Schrodinger equation (1.3). More 



precisely, the two x-parts are identical if Ux is identified with q, where u is the solution 



of (1.2) and q is the solution of (1.3). Although the corresponding t-parts are very 



different, this suggests a deep relationship between (1.2) and (1.3). Since all members 



of the hierarchy associated with an integrable equation admit the same x-part, it is 



reasonable to expect the equation describing the evolution of Ux according to ( 1.2 ) to 



be related to a member of the DNLS hierarchy. We will show that a simple change 



of variables indeed transforms (1.2) into the first negative member of this hierarchy. 



3.1. Gauge transformation. Our first step is to transform equation (1.2) by means 
of a gauge transformation. Replacing x by —x if necessary, we can assume that the 
parameters 7 and ly in (1.2) have the same sign. Then, letting a = > and 
b = 1/z/, the gauge tranformation 

\P\ 




^The variables A, z, and T in (|2.13|) are identified with u, t, and x in (|1.2|), respectively. The 
parameter v is identified with 



NONLINEAR PULSE PROPAGATION IN OPTICAL FIBERS 



7 



transforms (1.2) into 



utx — auxx = oib {—u + itylu] Ux), 
Introducing new variables by 

^ = X + ot, r = 

we arrive at the equation 



a = sgn p = ±1. 



(3.1) 



= u — ia\u\ u^, 



a = ±1, 



which is exactly (1.4) up to a relabeling of r and ^. Putting the various transforma- 
tions together we arrive at the following result. 

Proposition 3.1. Fix any nonzero values of the parameters z^, 7,p G M. Then 
u(x,t) satisfies equation (1.2) if and only if Ug{^,T) satisfies equation (3.1) with 
a = sgn^ujp) where u and Ug are related as 



t) = . /sgn(i^7p)^^e •^''u„{x-\ — t,- 



;t 



3.2. Bi-Hamiltonian structure. Equations (1.3) and (1.4) are equivalent to the 
systems 



(3.2) 

and 
(3.3) 



-zrr 



{q^r)x 



U — lUVUx 

V + iuvvx 



J1 



respectively, when r = aq and v = au. The operators 

dx\ J ^ / -qd-'^Q i + qdx^r\ 

, y ' ^ \-i + rd-^q -rd-^r ) ' 

form a compatible pair of Hamiltonian operators cf. [21] . Starting with the Hamil- 
tonian 



qrdx, 



the operators Ji and J2 generate a hierarchy of bi-Hamiltonian equations according 

tcEl 

.r , 



Kn = -^igradi?„„i = J2gradi?„, n G Z, 



where the infinite sequence of conservation laws {Hn}nez are constructed recursively 
from the relation 

Jigrad/^n = J2grad//ri+i, n G Z. 
The first few members of this hierarchy are presented in Figure [T] where we have 
written u and v for d~^q and d~^r, respectively. We immediately recognize equation 



The gradient of a functional F[q,r] is defined by 



grad F - 



6r ' 



wlienever there exist functions ^ and ^ such that, for any smooth functions ipi and ip2, 



d~e 



F[q + etpi,r + eip2] 



5F 3F . , 

-^<pi + -^'P2 ]ax. 
oq or 
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H2 



grad J {—qr) dx = gradi/o 




grad J (— nt> + ^uv'^Ux) dx = gradH-2 

Figure 1. Recursion scheme for the operators Ji and J2. The notation u = d~^q 
and V = d^'^r has been used when convenient. 



(3.2) as the second member and equation (3.3) as the first negative member of this 



hierarchy. 

Remark 3.2. Since the inverse of J2 is given exphcitly by 



'2 



-i + qd^ V qd^ ^q 



it is straightforward to implement the recursive scheme for construction of the H^s 
on a computer and thus obtain exphcit formulas for as many conserved quantities as 
one wishes. The H^s are local in q and r for n > 0, while they involve the nonlocal 
operator for n < 0. In both cases, the expressions grow rapidly in size, e.g. 
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and H4 are given by 

^3 = I { - ^Q'^r'^ - T^iq^r^r"^ + \qqxxr'^ + qqxTxr + 2q^rxxr + ^q^'rl + iqrxxx ]dx, 
Ha = j [ - ^q^r^ + ^qqlr^ - y i^'^r^^r^ + ^^^fer^ + hq^qxVxv'^ + y g^r^^r^ 

"-^q^rlr + \iqrxqxxr + ^iqqxTxxT + \iq^rxxxr + \iqqxrl 
b 2 2 2 2 

5iq Tx'^xx q^xxxx I dx. 



4. Traveling waves 



In this section we analyze the travehng-wave solutions of equation (1.2) taken in 
the simpler, but equivalent, form (1.4). Letting — > (at, fia:), we may assume 
that a = 1 in (1.4). Thus consider the equation 



(4.1) 

Substituting 



Utx = U — l\U\ Ux- 



u{x,t) = VF(x,t)e*^(^'*\ 
with real- valued functions W and x into (|4.ip, we obtain the coupled equations 

(4.2) 



-XxW^ - ixtXx + l)W + Wtx = 0, 
WxW^ + xtxW + Wtxx + XtWx = 0. 



We seek a solution of the form 

(4.3) x{x,t) = kx-nt + e{y), 



Wix,t) = ip{y), 



where k, J7, c are real parameters and y = x — ct. Equations (|4.2|) and (|4.3|) yield 
(4.4) 



{k + By) -{l-{k + Oy) {n + Cey)) if " Of yy = 



-Lfy [—if + ck + ^ + 2c6y) — c(p9. 



yy 







Multiply the second of these equations by ip and integrate the resulting equation to 
find 

{ck + n + 2c9y) Lf"^ + A = 0, 



1_ 

4 2 



for some constant of integration A. Solving for 9y^ we obtain 



(4.5) 



2ck(p^ - 2nif'^ + 4A 



We use this equation to eliminate 6y from the first equation in (4.4). Integration of 
the resulting equation multiplied by ^py yields 



(4.6) 



1 

16c2 
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Figure 2. The real part, imaginary part, and the absolute value at t = of the 
one-soliton solution u(x,t) given by equation The graphs correspond to the 

parameter-values A = 1, 7 = 1/2, and Sq = 0o = 0. 



where 



C3 

C2 
Cl 

CO 



ip^ + Csif^ + C2f^ + Ci(y9^ + Co 

4ck - 4n, 

4 (c^A;^ + n'^ + 2A + c(4 - 2kn)) , 

-32Bc, 

16^^ 



and B is anotlier integration constant. Solutions of (4.6) can only exist in regions 
where V{ip) < 0. Since cq > unless A = 0, we infer that solitary waves (i.e. waves 
such that (f decays to zero as ?/ — > ±00) can exist only when A = 0. A classification 
of the traveling waves can be obtained by analyzing the distribution of the zeros of V 
as the parameters c, k, A, B vary. Rather than completing this program in detail, 
we simply indicate how the one-soliton solutions arise and discuss the (non-)existence 
of peaked traveling waves. 



4.1. One-solitons. Equation (4.1) admits the one-soliton solutions [5] 



(4.7) 
where 
(4.8) 

(4.9) 



Us{x,t) 



A(e20+n + 1) 
t — 4A^a;) cos 7 



e 



2A2 

[t + 4xA'^) sin 7 



2A2 



+ eo, 
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and 

7 G (0, vr), A > 0, Sq G M, Qq G 
are four parameters, see Figure [2} In order to relate this family of one-solitons to 
the above description of traveling waves, we need to find parameters Cs,ks,^s, and 
functions 9s[y)-,^s[y) such that 

(4.10) Us{x, t) = ^s{x - c,t)e^('=^^-^^*+^»(^-^^*)\ 

The amplitude of Ug satisfies 



\Us 



2 sin^ 7 



A2(cos7 + cosh 29) ' 



This yields 

(4.11) vs{y) 



2 sin^ 7 



A2 (cos 7 + cosh (4A2y sin 7 + 20o)) ' 



4A4' 



It can be verified that the function (ps{y) given by (4.11 1 satisfies the ODE (4.6) 
provided that A = B = and 

cos 7 



(4.12) 
Substi 
(4.13) 



k, = 0, 



5^. 



A2 ■ 



Substitution of expression (4.11) for tps into ( |4.5[ ) gives after integration 
6siy) = — arctan ^tan tanh (2A^ysin7 + Qq 
-cot(7) (2A2ysin7 + eo) +0os, 



where Oqs is a constant of integration. We can determine 6*0^ by evaluating (4.10) at 
X = t = 0. This yields 



i\/2e<2^-*®«+^°+^''"*^''(*^'"(^)*^'"'"®°)+®°"°*^) Vcos7 + cosh(2eo) 
In summary, the one-soliton Ug can be written in the form (4.10) where Cs,k_ 



(4.14) 



6s{y),(ps{y) are given by equations (4.11 )-(4.14 ). This shows how the one-soliton 
solutions are related to the above traveling- wave analysis. 

4.2. Peakons. The Camassa-Holm equation is well-known for its peakons, which are 
solitons with a peak at their crest [14J (see Figure [s]). Since equation (1.2) is related 
to the NLS equation in a similar way that the Camassa-Holm equation is related to 
KdV, one may wonder if (1.2) also admits some kind of peaked solutio ns p| 

Let us first make some general comments on equations of the form (4.6) and the 



occurence of peakons. The ordinary differential equation (4.6) can be analyzed qual 



itatively by regarding (p{y) as describing the motion of a ball rolling in a potential 
well V{ip) with y representing time. Using this analogy, peaked solitary waves arise 
as follows. Consider a potential V{(p) with a qualitative shape as in Figure |4j Let 
the ball at time y = —00 be located at the crest of the hill at the point A. Give 
it a gentle push so that it starts rolling down the valley. If nothing intervenes, it 
will climb the opposite side of the valley until it reaches the point B where it turns 



■^More precisely, we say that a continuous real-valued function f{x) has a peak at xq if / is smooth 
locally on either side of xq and the left and right derivatives of / at xo are both finite but not equal. 
A peaked solution u{x, t) is a solution such that, for each fixed time t, u{-,t) has at least one peak. In 
the case of a complex-valued solution u{x, t), one could allow peaks in either the real or the imaginary 
part of u{x,t), or in both. 
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Figure 3. A peakon. The discussion in Subsection suggests that equation (1.2) 
admits no peaked solutions. 




Figure 4. An example of a potential well V{(p). 



around; it crosses the valley again and has just enough energy to climb all the way 
back up to the starting point A as y — > oo. This motion corresponds to a smooth 
solitary wave solution. However, if we intercept the solution before it reaches point B 
at some point C at which we let it (after a perfectly elastic collision) rebound, then 
it will climb up the first hill again and approach its starting position ^ as y ^ oo. 
The effect of the bounce is to create a peakon singularity at which the derivative (py 
jumps, t.py -ipy. 

Clearly, for most soliton equations there are plenty of ways to intercept solutions 
in this manner and thus obtain peaked solitons (see [22] for an example of this in 
the context of fiber optics). However, even though the peaked traveling waves so ob- 
tained satisfy the original PDE wherever they are smooth, the singularities in general 
introduce additional 5-functions so that the PDE is not satisfied in a reasonable weak 
sense. For most equations this rules out the existence of peaked solitons. On the 
other hand, if we consider the Camassa-Holm equation 



(4.15) 



Ut - Utxx + 3uUx = 2UxUxx + UUx 



the following happens. For a traveling wave u(x, t) 
yields 



ip{x — at), integration of (4.151 



(4.16) 



C(f + {C- (f)(pxx + 



1 2 « 

2^-- + 2' 



for some constant of integration a E M. In going from (4.15) to (4.16), the utxx and 
uUxxx terms have combined to give the term (c — ^)^xx- 

Moreover, (c - ip)ipxx is 



the only term in (4.16) which involves a second derivative of ^p. At a point where 
f has a peak, (pxx is proportional to a (5-function. Hence, the singularity-induced 
(5-function in (4.16 ) at a peak appears multiplied by c — ip, where c is the speed of the 
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wave. One could therefore hope for a peak to be allowed at a point xq provided that 
c — (p{xo) = 0, since this factor then kills the J-function. This is indeed what happens: 
a peak in a traveling wave of the Camassa-Holm equation is admissible exactly when 
the height of the peak equals the speed of the wave, i.e. when ip = c at the crest cf. 



On the contrary, we will now argue that for equation ( 1.2 ) the (5-functions produced 



by the x-derivatives of a peakon candidate stand alone and cannot be cancelled by 



other factors. For simplicity we consider the equivalent, but simpler, equation (4.1); 



the argument for (1.2) is analogous. When substituting a traveling-wave solution of 



the form (1.5) into equation (4.1), the terms involving a second derivative oi ip oi 9 
can be read off from equation (4.4). They are 



—cipyy in the equation for the real part, 
—apOyy in the equation for the imaginary part. 

For the solution to have a peak at yo at least one of the functions (pyy or 6yy must 
involve a (^-function at yo. The only way for the real and the imaginary parts to 
vanish simultaneously is therefore that either (i) c = or (ii) only the function 9 
has a peak at yo and ip{yo) = fyyivo) = 0. The following observations indicate that 
neither of these two options is viable. 
Regarding (i) we note that if c 



then the second equation in (4.4) becomes 

-<fy{-ip'^ + n) = 0, yeR. 

By the assumption that the amplitude of a peakon is a continuous function, this 
shows that the amplitude (p is in fact constant. It then follows from the first equation 
in (4.4) that 9y is also a constant function given by 

9y = 



Thus, no peakons arise when c = 0. 

Regarding (ii) we note that a vanishing amplitude(^(yo) = at yo implies that the 
phase 9(yo) is undetermined a priori. Of course 9 may contain peaks (or jumps for 
that matter) in any intervals where ip = 0, but these peaks are unobservable and are 
just a consequence of the breakdown of the polar coordinate system at the origin. 



As soon as (p becomes nonzero the value of 9y is fixed by the equations in (4.4) and 
cannot jump. Hence, neither of the options (i) or (ii) leads to existence of peakons. 
The above discussion indicates that (4.1) admits no peaked solutions within the 



class of traveling waves of the form ( |l.5[ ) . Although one could imagine that there exist 
peaked solutions of a more general type, similar arguments as above applied directly to 
equation ( |4.1[ ) make this unlikely (since utx is the only term in the equation involving 
two derivatives, there is no other term to balance the 5- function generated by utx at 
the peak). In conclusion, it appears that equation (1.2) admits no peakons which are 
weak solutions in a reasonable sense. 



References 

[1] G. P. Agrawal, Nonlinear Fiber Optics, fourth edition, Academic Press, 2007. 

[2] V. E. Zakharov and A. B. Shabat, Exact theory of two-dimensional self-focusing and one- 
dimensional self- modulation of waves in nonlinear media, Soviet Phys. JETP 34 (1972), 
62-69. 

[3] A. S. Fokas, On a class of physically important integrable equations, Phys. D 87 (1995), 
145-150. 



14 



JONATAN LENELLS 



[4] B. Fuchssteiner and A. S. Fokas, Symplectic structures, their Backlund transformation 

and hereditary symmetries, Phys. D 4 (1981), 47-66. 
[5] J. Lenells and A. S. Fokas, On a novel integrable generahzation of the nonhnear Scrodinger 

equation, submitted. 

[6] F. Calogero, Why are certain nonhnear PDEs both widely apphcable and integrable?, in: 

What is integrability? (V.E.Zakharov, editor). Springer, 1990, 1-62. 
[7] Y. Kodama, Optical solitons in a monomode fiber, J. Stat. Phys. 39 (1985), 597-614. 
[8] Y. Kodama and A. Hasewaga, Nonlinear pulse propagation in a monomode dielectric 

guide, IEEE J. Quantum Electron. 23 (1987), 510-524. 
[9] D. Anderson and M. Lisak, Nonlinear asymmetric self-phase modulation and self- 
steepening of pulses in long optical waveguides, Phys. Rev. A 27 (1983), 1393-1398. 

[10] N. Tzoar and M. Jain, Self-phase modulation in long-geometry optical waveguides, Phys. 
Rev. A 23 (1981), 1266-1270. 

[11] E. Mjolhus, On the modulational instability of hydromagnetic waves parallel to the 
magnetic field, /. Plasma Phys. 16 (1976), 321-334. 

[12] D. J. Kaup and A. C. Newell, An exact solution for a derivative nonlinear Schrodinger 
equation, /. Math. Phys. 19 (1978), 789-801. 

[13] H. P. McKean, The Liouville correspondence between the Korteweg-de Vries and the 
Camassa-Holm hierarchies. Comm. Pure Appl. Math. 56 (2003), 998-1015. 

[14] R. Camassa and D. D. Holm, An integrable shallow water equation with peaked solitons, 
Phys. Rev. Lett. 71 (1993), 1661-1664. 

[15] J. Hermmann and B. Wilhelmi, Lasers for ultrashort light pulses, North- Holland, 1987. 

[16] H. H. Chen, Y. C. Lee and C. S. Liu, Integrability of nonlinear Hamiltonian systems by 
inverse scattering method, Phys. Scr. 20 (1979), 490-492. 

[17] S. Yu. Sakovich, Integrability of the higher-order nonlinear Schrodinger equation revis- 
ited, arXiv:so lv-int/990 6012v l. 

[18] C. Gilson, J. Hietarinta, J. Nimmo, and Y. Ohta, Sasa-Satsuma higher order nonlinear 
Schrodinger equation and its bilinearization and multi-soliton solutions, Phys. Rev. E 68 
(2003), 016614. 

[19] R. Hirota, Exact envelope-soliton solutions of a nonlinear wave equation, J. Math. Phys. 
14 (1973), 805-809. 

[20] N. Sasa and J. Satsuma, New-type soliton solution for a higher-order nonlinear Schrdinger 

equation, J. Phys. Soc. Jpn. 60 (1991), 409-417. 
[21] W. Ma and R. Zhou, On inverse recursion operator and tri-Hamiltonian formulation for 

a Kaup-Newell system of DNLS equations, J. Phys. A 32 L239-L242. 
[22] L. Cheng, The cusp sohtons in optical fibres, J. Phys. A, 29 (1996), 141-142. 
[23] J. Lenells, Traveling wave solutions of the Camassa-Holm equation, J. Diff. Eq., 217 

(2005), 393-430. 

J.L.: Department of Applied Mathematics and Theoretical Physics, University of 
Cambridge, Cambridge CBS OWA, UK 

E-mail address: j.lenells@damtp.cam.ac.uk 



